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Q\ , Abstract 
O 



The skewmg factor, defined as the ratio of the imaginary parts of 
the amplitudes ImA{'y* p — > 'y* p) / Im Ai^-y* p jp) is extracted for 
■ the first time from the data using recent Deeply Virtual Compton 

Scattering (DVCS) and the inclusive inelastic cross section measure- 
ments at DESY-HERA. The results values are compared to theoretical 
predictions for NLO QCD and the colour dipole approach. 



Introduction 

X 

5t ! It is well known that the cross section of hard scattering processes can be 

described as the convolution of parton distributions (PDFs) and the cross 
sections of hard subprocesses computed at parton level using perturbative 
QCD. The usual PDFs, obtained from inclusive experimental data, is the 
diagonal element of an operator in the Wilson operator product expansion 
(OPE). On the other hand, there is a set of exclusive reactions which are 
described by the off-diagonal elements of the density matrix, where the mo- 
mentum, helicity or charge of the outgoing target are not the same as those 
of the corresponding incident particle. Examples of such reactions are the 
virtual photon Compton scattering (DVCS) (7*^ '-fp) and the vector me- 
son electroproduction (7*^ — > Vp). In these cases, the difference with the 
inclusive case is the longitudinal components of the incoming and outgoing 
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Figure 1: Amplitude level diagrams for DIS (a) and DVCS (b) 

proton momentum, which depend on the photon virtuality and the 'j*p 
center of mass energy W. 

Precision data are becoming available for hard scattering processes whose 
description requires knowledge of these off-diagonal (or skewed) parton distri- 
butions. The exclusive diffractive DVCS process at large provides a com- 
paratively clean procedure for extracting information on the gluons within 
the proton in a non-forward kinematic case. In recent years the study of Deep 
Inelastic lepton-proton Scattering (DIS) has produced detailed information 
of the proton structure in terms of Parton Distribution Functions (PDFs). 
The optical theorem states that the leading DIS process of a single virtual 
photon exchange can be viewed as identical to the forward elastic scattering 
of a virtual photon from the proton. This scattering, known as the forward 
Compton scattering process, is the same as a DIS process (see Fig. HJ-a) with 
its mirrored reaction and relates the total cross section atot{l*P X) with 
the imaginary part of the forward amplitude ImA{'~f*p 7*p). 

The forward Compton scattering 'y*{q)p{p) l*{<l)p{p) can be general- 
ized to the non forward case (at non-zero angles): 'J*{q)p{p) l^*\(l2) p{P2) 
with q2 ^ Q and P2 ^ P and being experimentally accessible through the so 
called Deeply Virtual Compton Scattering (DVCS) process, where a real pho- 
ton is observed in the final state (see Fig. ^b). In analogy with DIS which 
relates the imaginary part of the forward Compton amplitude and PDFs, 
in the DVCS process parton distributions are also used to parameterise the 
non-perturbative part of the interaction in the non-forward region. In this 
last case, the distributions have to be calculated in terms of non-diagonal, 
rather than diagonal, PDFs. They are commonly called Generalised Parton 
Distributions (GPDs). They reduce to the ordinary PDFs in the limit of zero 
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skewedness, i.e., the forward limit. 

Let us first summarize the main theoretical approaches describing the 
DVCS process. From the perturbative QCD point of view, GPDs have been 
studied extensively and incorporate the partonic and distributional behavior. 
They are two-parton correlation functions that allow accessing parton corre- 
lations within hadrons. For fixed skewedness, they are continuous functions 
and cover two distinct regions, the DGLAP and the ERBL regions, in which 
their evolution with the hard scale is driven by the corresponding evolution 
equations. Moreover, the Lorentz structure of the GPDs implies a polyno- 
miality condition; namely, their (n — l)-th moments are polynomials in the 
square of the skewedness of degree smaller than n/2. The most promising 
approach for GPD parameterisations relies on setting them equal to the for- 
ward PDFs using a forward model with suitably symmetrized input GPDs 
in the ERBL region constructed in order to satisfy polynomiality for the first 
two momenta |7I . These theoretical studies reproduce satisfactorily the main 
features of DVCS data jj]. On the other hand, a simple and intuitive high 
energy formalism has recently been used to describe DVCS jH]. It is given by 
a colour dipole picture, which allows extending in the studies of saturation 
physics at low virtualities of the incoming photon. This is an advantage in 
comparison with the pQCD formalism which is limited by the initial scale for 
QCD evolution to the order of Qo ~ 1 GeV. Additionally, the roles played 
by QCD evolution and skewing effects have been recently shown in the scope 
of this formalism. The off-diagonal effect has been investigated through a 
simple phenomenological parameterisation. The overall picture is in good 
agreement with the experimental data ^llj. 

In this work, we will address the question of the extraction of the skewing 
correction using the available data on DIS and DVCS. This can be performed 
by computing the ratio of the imaginary parts of the amplitudes for these 
reactions. The skewing factor R will depend on the experimental measure- 
ments and therefore possess a small degree of model dependence. We there- 
fore have the attractive possibility in determining its dependence upon 
virtuality and upon W energy. This knowledge will enable us to gain insight 
into the skewedness correction to the non-forward observables. Namely, the 
non-forward cross section and distributions can be easily obtained from the 



3 



forward ones by simply multiplying them by the skewing factor. This has 
direct implications in computing the exclusive meson and heavy boson pro- 
duction via the forward scattering amplitude, which is better constrained 
from the inclusive measurements. 

The paper is organized as follows; in the next section, the skewing factor 
R is defined and extracted from experimental DESY-HERA results on DIS 
and DVCS. In section |21 we perform a comparison of the skewing factor with 
two different theoretical approaches, the perturbative NLO QCD calculation 
and the colour dipole formalism. A qualitative analysis of the extracted 
skewing factor is then performed in a colour dipole approach. In the last 
section we summarize the main results and discussions. 

1 Definition and extraction of the skewing 
factor R 

Lets define a basic quantity giving an overall measurement of the skewing 
properties, which includes both the non-forward kinematics and the non- 
diagonal effects. Namely, we set the ratio between the imaginary parts of the 
DIS and DVCS (forward) scattering amplitudes at zero momentum transfer: 





^ 7* + 


t=o 


Im A (7* + p - 


7 +p) 





(1) 



where t is the square of the four-momentum exchanged at the proton vertex. 

The scattering amplitude for the DIS process can be directly obtained 
from the DIS cross section and experimentally measured at DESY-HERA, 
that is ImA{YP ~^ 1*P) ~ o'tot{'y*P — > X)- In fact, the DIS amplitude 
can be written down using the usual pQCD fits for the proton structure 
function, o"7ot^ = {4:n'^a/Q^) F2{x,Q'^). The DVCS scattering amplitude can 
be obtained from the recent measurements on the Deeply Virtual Compton 
Scattering cross section. In this case, the t dependence of the amplitude 
can be assumed to be factorised out and parameterised as an exponential. 
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oc exp(— and the total DVCS cross section can be related to the corre- 
sponding amplitude at t = 0; namely, 



t=o 



(2) 



167r6 

where h is the t slope parameter referred to above. The expression above 
can be corrected by taking into account the contribution from the real part 
of the amplitude by multiplying Eq. (j2I) by a factor (1 + rf), where 7] is 
the ratio of the real to imaginary part of the DVCS amplitude. The typical 
contribution for the kinematic window available at HERA is of the order of 
20 %. In particular, to a good approximation, r] can be calculated using 
dispersion relations to ~ tan(7rA/2), where A = A((5^) is the effective 
power of the Bjorken x dependence of the imaginary part of the amplitude. 
For an estimate of the real part contribution we use the effective power for 
the inclusive deep inelastic reaction taken from Ref. jS]. 

Considering the calculation discussed above, we can rewrite the skewing 
factor as a function of the total cross sections for DIS and DVCS. The factor 
reads as, 

a(7*p^X) ^{l+n^) v^aF|'(x,g2)^(l + r^2) 
^ / = ~~ry2 / • l<jj 

Main theoretical uncertainties come from the b slope and from the estimate 
of the real part contribution. In our further calculation, one uses the b value 
extracted from the recent measurements of the t-dependence of DVCS data. 



To extract the factor R factor, Eq. Q, we use recent DVCS measure- 
ments at HERA PI 12] and the DIS cross cross section is obtained using the 
pQCD fits of the Fj* structure function from Ref. 1^. The factor R is shown 
as a function of in Fig. |21 and as a function of energy W in Fig. El The 
corresponding values are given in tables ^ and El The inner error bars rep- 
resent the statistical error. The full error bars is the quadrature sum of the 
statistical, systematic and normalisation (precision of the b measurement of 
HI) uncertainties^. A dependence is observed, decreasing from i? ~ 0.7 

^important contribution to systematic effects in the HI measurement are not considered 
in the ZEUS measurement (e.g. the error on the Bethe-Heitler subtraction whicli dominates 
the fuU systematic error at large W values) 
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Figure 2: The skewing factor R as a function of Q^atW = 82 GeV. The 
points correspond to the R extraction apphed to DVCS measurements of 
HI (bullets) and ZEUS (triangles). The curves represent R extracted from 
Freund et al. prediction based on MRST 2001 and CTEQ6 PDFs. 



for Q2 ^ 2 GeV2 down to R - 0.3 for ~ 85 GeV^. Qualitatively, the data 
seem to be consistent with a dependence like R oc 1/logQ^. An almost flat 
W dependence is observed within the present precision. This feature can be 
easily understood by inspecting Eq. (jHl), since the W dependence of both 
the DIS and DVCS cross section is power-like having a proportional effective 
power. Namely, a^is oc VT^^ and ctdvcs W^^. The mean value R ~ 0.5 is 
consistent with its early theoretical estimates using the aligned jet model jH] 
and the colour dipole saturation model 0. 



Q' [GeV2] 


R using HI data 




[GeV^] 


R using ZEUS data 


3 


0.62 ±0.05 ±0.07 




7.5 


0.53 


±0.02 ±0.03 


5.25 


0.66 ±0.06 ±0.09 




12.5 


0.48 


±0.02 ±0.02 


8.75 


0.57 ±0.04 ±0.07 




20 


0.45 


±0.03 ±0.02 


15.5 


0.55 ±0.05 ±0.08 




32.5 


0.40 


±0.05 ±0.03 


25 


0.45 ±0.06 ±0.06 




55 


0.41 


±0.09 ±0.05 


55 


0.42 ±0.07 ±0.07 




85 


0.28 


±0.09 ±0.04 



Table 1: The skewing factor R for different values at W = 82 GeV 
according to the DVCS cross section measurement of HI (left table) and 
ZEUS (right table). 
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Figure 3: The skewing factor R as a function of W at Q"^ = 8 GeV^. (see 
caption of Fig. 0). 



W [GeV] 


R using HI data 


45 


0.56 ±0.03 ±0.05 


70 


0.60 ±0.04 ±0.06 


90 


0.59 ±0.04 ±0.07 


110 


0.60 ±0.07 ±0.10 


130 


0.53 ±0.10 ±0.11 



W [GeV] 


R using ZEUS data 


45 


0.50 


±0.03 


±0.02 


55 


0.47 


±0.02 


±0.02 


65 


0.46 


±0.02 


±0.02 


75 


0.47 


±0.03 


±0.02 


85 


0.52 


±0.03 


±0.02 


95 


0.57 


±0.05 


±0.03 


105 


0.43 


±0.03 


±0.02 


115 


0.51 


±0.04 


±0.02 


125 


0.53 


±0.05 


±0.02 


135 


0.51 


±0.04 


±0.02 



Table 2: The skewing factor R for different W values at = 8 GeV^ 
according to the DVCS cross section measurement of HI (left table) and 
ZEUS (right table). 
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2 Theoretical Predictions 



In this section we compare the results extracted from experimental data to 
different theoretical predictions. First, we contrast them with the perturba- 
tive QCD approach at NLO accuracy. This prediction is dependent on the 
Generalised Parton Distributions (GPD), which have been parameterised in 
Ref. [7j and applied to describe the recent DVCS data OEI- The second 
theoretical approach [HI E] is given by the colour dipole approach, where the 
basic building blocks are the photon wavefunctions and the dipole cross sec- 
tion. The non-forward kinematics is encoded by the wavefunctions whereas 
the off-diagonal effects can be built-in in the parameterisations for the dipole- 
target cross section. 

2.1 NLO QCD Predictions 

The DVCS cross section has been calculated at NLO in perturbative QCD by 
Freund and McDermott [21 El using two different GPD parameterisations [7] . 
The MRST2001 and CTEQ6 parameterisations of the PDFs are used in the 
DGLAP region and polynomial form is used in the ERBL region, ensuring a 
smooth continuation between the two regions. Both the skewing (^) and the 
dependence are generated dynamically. The t dependence is factorised 
out and assumed to be e"**'*'. In its original form, the approach above used 
a Q^-dependent h slope. We have verified that the recent data on DVCS [T] 
can be reasonably described using a fixed h slope. 

Applying the same method presented here (i.e. using the same for 
the inclusive cross section, the values obtained are shown as a function of 

in Fig. 121 and as a function of W in Fig. |21 The value of 6 = 6.02 ± 
0.35 ± 0.39 GeV~^ measured by HI [1^ has been used in all predictions. The 
error on the h value is not applied to theoretical predictions as it is already 
included in the total error of the data points. A good agreement between 
these predictions and the data points is found, describing well the absolute 
value of R and its kinematic dependences. The NLO calculation presents no 
W dependence. The result normalisation is strongly dependent on the choice 
for the gluon distribution and the deviations seem to be amplified at lower 
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Q^. Concerning the behavior on Q^, it seems that for values below ~ 4 
GeV^ the DGLAP evolution starting at = 1 GeV^ has too little phase 
space to fully generate the gluon distributions (dotted part of the curves in 
Fig. 121). The underestimation of R for virtualities below 4 GeV^ may be due 
to a larger value for the DVCS cross section calculated at NLO where the 
gluon distribution is underestimated as an effect of the too small phase space 
for evolution. 

2.2 Colour Dipole Model Predictions 

In the colour dipole approach, the DIS (or DVCS) process can be seen as a 
succession in time of three factorisable subprocesses: i) the virtual photon 
fluctuates in a quark-antiquark pair, ii) this colour dipole interacts with the 
proton target, iii) the quark pair annihilates into a virtual (or real) photon. 
The imaginary part of the DIS (or DVCS) amplitude at t = is expressed 
in the simple way P, 

1 oo 

TmA{W,Qi,Q2) = J2 [ ^zf c/V^;._^(z, r, Qj) adip{x,r)<i/T,L{z, Qj) , (4) 

where \l/(2, r, (^1,2) are the light cone photon wave functions for transverse 
and longitudinal photons. The quantity Qi is the virtuality of the incoming 
photon, whereas Q2 is the virtuality of the outgoing photon. In the DIS case, 
one has = = and for DVCS, Qj = and Ql = 0. The relative 
transverse quark pair (dipole) separation is labeled by r whilst z (respec. 
^ — z) labels the quark (antiquark) longitudinal momentum fraction. 

It should be noticed that the non-forward kinematics for DVCS is en- 
coded in the colour dipole approach through the different weight coming 
from the photon wavefunctions in Eq. (jl)). The off-diagonal effects, which 
affect the gluon and quark distributions in the pQCD approaches, should be 
included in the parameterisation of the dipole cross section. At the present 
stage of the development of the colour dipole formalism, we have no accu- 
rate theoretical arguments on how to compute skewedness effects from first 
principles. A consistent approach would be to compute the scattering am- 
plitude in the non-forward case, since the non-forward photon wave function 
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has been recently obtained in Ref. [^. In this case, the dipole cross section, 
o'dip{xi,X2,r,A), depends on the hght cone momenta xi and X2 carried by 
the exchanged gluons, respectively, and on the total transverse momentum 
transfer A. In this case, additional information about the dependence upon 
A is needed for the QCD Pomeron and proton impact factor. The forward 
dipole cross section is recovered at xi = X2 and A = 0. 



In Ref. jTT] an estimate of the skewedness for the dipole cross section 
has been performed using the approximation of Ref. ^3] where the ratios of 
off-diagonal to diagonal parton distributions are computed. The behavior of 
those ratios are given explicitly by: 

where p = for quarks and p = 1 for gluons, whereas A is the effective 
exponent of the parton distribution. The ratio is larger for singlet quarks 
than for gluons. In the colour dipole picture, the DVCS cross section is driven 
mainly by gluonic exchanges. In our numerical computations, we use A = 
XiQ"^) as obtained from the DVCS scattering amplitude and the skewedness 
effect is given by multiplying the original (without skewedness) total cross 
section by the factor RglQ"^). A different implementation of the skewedness 
correction is to use the approximation x = 0.41 xbj, which produces the same 
effects as for Rg at the presently covered kinematic range and precision of 
the measurement. This different procedure is particularly suitable for dipole 
cross sections which are parameterised by using a gluon distribution, namely 
adip oc r^xg {x,c/r^). 

Using this colour dipole prediction (with DGLAP evolution and off-diagonal 
skewing factor Rg), one can again extract the R factor, as shown as a func- 
tion of in Fig. |3] and as a function of W in Fig. A good agreement 
with points extracted from the data is found. The curve obtained without 
applying the off-diagonal skewing factor Rg is also presented, enabling one to 
isolate the two effects: the non-forward kinematic and the off-diagonal gluon 
distributions. En passant these results demonstrate that the approximation 
of Eq.Elis reliable within the present precision. Furthermore, the off-diagonal 
effects contribute mostly for the overall normalisation of the skewing factor 
whereas the behavior on seems to be driven by the non-forward kinematic 
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Figure 4: The skewing factor R as a function of Q^atW = 82 GeV. The 
curves represent the theoretical predictions of the colour dipole approach (see 
text) with (solid) and without (dashed) applying the off-diagonal skewing 
factor Rg. 



effect. Tfiis prediction exliibits a moderate increasing of R witli W wliicli is 
not present in tlie NLO prediction. 

Let us perform a qualitative analysis in order to obtain a pliysical picture 
for tfie W and dependences of the skewing factor. We take the colour 
dipole approach and the saturation model for the dipole cross section [Tn] . 
considering for simplicity the massless case. The final result should be suffi- 
ciently independent of this particular choice for the dipole cross section. One 
can use an approximated form for the dipole cross section from the satura- 
tion model; adip = a^r"^ Ql{x) / A (colour transparency) for < A/Ql{x) and 
o"dip = ctq (black disk hmit) for > A/Ql^x). The saturation scale is given by 
Ql{x) oc Aqcds;"'^ and gives the onset of nonlinear corrections to the QCD 
dynamics. As the DVCS data are predominantly at intermediate Q^, we will 
be interested in small size colour dipoles (here, the dipole cross section is 
dominated by colour transparency behavior). In this particular region, we 
found in Ref. [Hj the qualitative behavior for the total DVCS cross section, 
taking into account only the non-forward kinematic effect. 



Q4 

f^(7> ^ 7P) oc ^ 




(6) 
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Figure 5: The skewing factor R as a function of W at Q'^ = 8 GeV^. (see 
caption of Fig. Ej). 

where the large logarithm in Eq. © originates from the intermediate region 
2/Q < r < 2/Qs{x) and the remaining comes from the region 2/Qs{x) < r < 
2/Q m the dipole size integration in Eq. 

The qualitative result above can be contrasted with the DIS case, where 
the 7*p total cross section grows as ~ W'^'^/Q'^ at intermediate Q^. In par- 
ticular, DIS data are proven to exhibit a geometric scaling pattern in the 
scaling variable r = Q'^/Ql- It is also easy to demonstrate that DVCS 
data also should present such a scaling behavior. In particular, at inter- 
mediate it is shown that the experimental DIS cross section behaves as 
a{'j*p X) ~ 1/r = Ql/Q"^. Putting this information together and the 
result for the DVCS cross section presented in Eq. (jH)) into the definition of 
the skewing factor, Eq. (jH)), one obtains. 



where in the DESY-HERA domain the saturation scale takes values of = 
1 — 2 GeV^. Therefore, the qualitative behavior for the skewing factor 
takes the form R oc 1/(1 + logQ^). On the other hand, one has ~ 
x^^ {W/QY^ which at fixed virtualities the skewing factor behaves as oc 
1/(1 — e log W). These qualitative results are in agreement with the behavior 
of the experimental extraction of R. One notice that the early theoretical 
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determinations for R predict R ~ log"^[l + (Q'^/M^)], with ~ 0.4-0.6 
GeV^, in close proximity with the present estimate [S]. 



3 Summary 



The skewing factor R, defined as the ration of the imaginary parts of the 
DIS and DVCS amphtudes has been extracted from experimental data for 
the first time. In its determination one uses the recent Deeply Virtual Comp- 
ton Scattering (DVCS) and the inclusive inelastic cross section measurements 
at DESY-HERA. The main theoretical uncertainties come from the b slope 
and from the estimate for the real part contribution. One founds a de- 
pendence qualitatively consistent with the form R cx 1/logQ^. No W de- 
pendence is observed within the current precision. The mean value R ~ 0.5 
is consistent with previous theoretical estimates in the jet aligned model and 
colour dipole picture. The extracted R is contrasted to theoretical predic- 
tions of NLO QCD and colour dipole approaches. The first one uses two 
different GPD parameterisations and produces a consistent description for 
the W and dependences of the skewing factor. The colour dipole formal- 
ism, supplemented by a dipole cross section including QCD evolution and 
off-diagonal corrections successfully describes R. In particular, a qualitative 
understanding of the dependences of R is obtained relying on simple argu- 
ments in the dipole formalism. This study confirms a behavior of the form 
R oc 1/(1 + logQ^) and flat dependence upon W. Beyond the successfully 
described dependences in and W achieved here, to perform a further dis- 
cussion on the normalisation of R, a better agreement within experimental 
DVCS cross section measurements has to be achieved. 
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